Detection and typicality of bound entangled states 
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We derive an explicit analytic estimate for the entanglement of a large class of bipartite quantum 
states which extends into bound entanglement regions. This is done by using an efficiently com- 
putable concurrence lower bound, which is further employed to numerically construct a volume of 
3x3 bound entangled states. 
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It is one of the most challenging and fundamental is- 
sues in quantum information science to decide whether 
a given quantum state can exhibit quantum correla- 
tions, i.e., whether it is entangled. This question is 
fundamental inasmuch as it rephrases the quest for the 
quantum-classical demarcation line, and it is also of po- 
tentially enormous practical relevance ~ in view of the 
many applications of quantum theory in modern infor- 
mation technology. Since entanglement is fragile, hard 
to screen against the detrimental influence of decoher- 
ence and rapidly reduced to a residual level under envi- 
ronment coupling, it is important to realize that even 
quantum states which are "close" to separable states 
and, in this sense, carry only residual amounts of en- 
tanglement, still might be used to accomplish typical 
tasks of quantum information processing, after "distilla- 
tion" pLj : Many weakly entangled states can be processed 
to condense their collective entanglement content in one 
strongly entangled state, which then can be used to solve 
the predefined task. 

However, there are entangled states from which no en- 
tanglement can be distilled, accordingly called bound en- 
tangled states [2|. If, e.g., two parties were to set up a 
quantum channel by sharing an entangled state, environ- 
mental noise can escort that state to a bound entangled 
state, thus preventing any subsequent distillation - the 
quantum communication channel will be ill-fated. An 
entangled state p that is positive under partial transpose 
(PPT) - i.e. (l(8iT')(py_fei) = Pii.kj > - was shown to be 
bound entangled j3J|U[S]. In other words, the only eas- 
ily computable entanglement measure to date [6j fails to 
detect exactly those states that represent a severe prob- 
lem for quantum communication protocols. It is there- 
fore mandatory to develop tools to efficiently map out 
the volume of bound entangled states, which is hitherto 
barely characterized: So far, only continuous families of 
optimal entanglement witnesses could be used to delimit 
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bound entangled states [Tj [3 El UHl HI] , by intersecting 
the volume of entangled states detected by the witnesses 
with the volume of quantum states with positive partial 
transpose. All quantum states within this intersection 
are bound entangled. The crux of this method lies in 
the difficulty of constructing optimal witnesses, which is 
known to be a computationally hard task, in general. 
Here we show that an algebraic lower bound of entan- 
glement when quantified by concurrence - which can be 
evaluated analytically or by numerical diagonalization - 
can be employed for efficient detection of an important 
class of bound entangled states of finite dimensional bi- 
partite quantum systems. 

Let us start by a short recollection of the basic defi- 
nitions of concurrence and its lower bound as employed 
hereafter. As shown elsewhere [12 [T^], Wootter's origi- 
nal concurrence definition for pure states |14j can be re- 
expressed in terms of the expectation value of a projector- 
valued operator A, 

C{-^) = , (1) 

where A acts simultaneously onto two versions of the 
state: 

^ = 4 ^ {\ik3l)-\3kil)-\iljk) + \jlik))-{h.c.) (2) 

i<j,k<l 

Here, i and j enumerate the basis vectors of the first par- 
tition, and k and I the second partition's. This definition 
of concurrence can be generalized for mixed states p, as 
an infimum over all possible pure state decompositions 
defined in terms of probabilities pi and pure states 

CiP) ^ , inf Vp,(vE',|®(vE-,|^|vI/,)0|vI/,)V2.(3) 

This latter optimization problem has an explicit algebraic 
solution for pairs of qubits [15' , but admits only numer- 
ical solutions or algebraic estimates if the system size 
is increased - either by an increase of the constituents' 
sub-dimension, or of their number. 
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In the following, we will use a specific, algebraic esti- 
mate, the quasi pure lower bound, which is easily evalu- 
ated (by diagonalization of a matrix of the same dimen- 
sion as p) and known to yield good estimates for weakly 
mixed states [16 . In short, it is obtained from the sin- 
gular values Si of a matrix with elements 

%j=^JIuI]{^,\<»{^j\x), (4) 

that can easily be constructed with the spectral de- 
composition p = X^i ™d choosing \x) oc 
y^l^'o) (8i 1^*0) , with I'^o) the dominant eigenvector of p 
(associated with the density matrix' largest eigenvalue). 
The concurrence can then be bounded from below by 



C (p) > Cqp (p) max 




for arbitrary states p (in contrast to witnesses, which 
need to be tailored for the detection of specific states). 
Therein, Sq denotes the largest singular value of ma- 
trix T. We will employ this quasi pure lower bound Cqp 
throughout the sequel of this paper. 

We now set out for identifying a volume of bound 
entangled states within the set of d dimensional Bell- 
diagonal states, a class of states of special importance 
for quantum key distribution protocols [17], and whose 
structure is not completely known [TH]. These are given 
as convex sums of maximally entangled Bell-like states, 

d-l 

P = ^ Afc/Pfci , (6) 

fc,/=0 



l2(d-l) 



which can be plugged into (|5| to obtain the desired re- 
sult (note that the singular values now carry four indices: 
the two upper-indices refer to the Bell state with largest 
eigenvalue, and the two lower indices are the labels of the 
Bell basis) . This represents the first analytical estimation 
of concurrence for a family of states that encompasses 
bound entangled states. 

We now apply this result to delineate the area of bound 
entangled states within the class of 3 x 3 Bell diagonal 
"line states" defined as (see also [51 [20]! 

P = ^ ~ " ~ ^ ~ ^ 1 + 7 -Poo + a ^10 + /3 ^20 ■ (9) 

(The same Weyl operator generates \^lio) from |rioo), as 
\il2o) from |riio) - therefore these states extend along a 
"line" .) For these states, the existence of a non- vanishing 
area of parameter space giving rise to bound entangled 



with probabilities Xki > 0, "^ki^ki = 1, and P^i = 
\^ki){^ki\ the projectors onto the Bell states 

l"«) = ;7^E^'^''l^)l-^ + 0- (7) 

The latter are transformed into each other by local uni- 
tary operations, e. g. by the Weyl operators Wki — 
+ such that \ flki) = (Wki^tmoo)- 

From several copies of the Bell state |iloo), a Bell- 
diagonal state is generated by introducing simple errors 
such as phase-shifts and bit-translations, and error cor- 
rection in general allows to reverse this process, and thus 
to distill a maximally entangled Bell state from a "reser- 
voir" of Bell diagonal states. Nonetheless, bound entan- 
gled Bell-diagonal states which do not admit entangle- 
ment distillation do exist [T5], and, as we will show in 
the following, can be effortlessly detected. 

For this purpose, we first derive Cqp for arbitrary Bell 
diagonal states as defined in ([6|. To do so, suppose that 
Pnm has the largest weight in ([6]), i.e., that p be quasi 
pure with respect to Pnm- We then can construct the ma- 
trix T in Q with the choice |x) (x A\^nm) \^nm), and 
the singular values of T are given by the square roots of 
the eigenvalues of TT^ - which itself can be shown to be 
a Bell diagonal matrix. Consequently, the singular values 
can be readily read off from 5^;" = {flki\%imX!m\^kiV^^ , 
with the explicit expression 



states had already been demonstrated through the op- 
timization of witness operators [HI UHl HI]- With the 
present approach, we can effortless scan the entire a-/3 
plane for fixed 7. The results, for 7 = in Fig. [T] show 
that the intersection of the area of positive Cqp with the 
area of positive partial transpose perfectly reproduces the 
area identified by the witness approach. Also note that 
the quasi pure bound here provides fully reliable informa- 
tion, despite the fact that the identified bound entangled 
states are rather mixed (with purities somewhere around 
0.17, and a minimum value 1/9 ~ 0.111). Fig. |2]shows 
analogous results for 7 < 0. 

As a byproduct of these results, we could exactly 
parametrize the borderline of bound entangled states and 
of those with positive partial transpose via the witness 
approach W. For positive a (which suffices due to the 
apparent symmetry of parameter space as spelled out by 



V2n-fc)(mod d),(2m-;)(mod d) 



1/2 



(8) 
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FIG. 1: (Color online) (a) Parameter space of the "line" states, as defined in for 7 = 0: The crosshatched dark gray area 
depicts states with positive partial transpose, while the quasi pure lower bound Cqp, Eq. (|5|, vanishes in the light gray area. The 
black area corresponds to states that have positive partial transpose and Cqp exhibit positive values, and thus is associated with 
bound entangled states. The solid gray triangle encompasses the positivity region (p > 0). The curves refer to the borderlines 
in Eq. ( 10 1: The full line delimits the positive partial transposition region, and the dashed ellipses the bound entangled region, 
when over the PPT region. Note that some entangled states, with negative partial transposition, have Cqp = 0, but this was 
not observed in areas of bound entanglement, (b) Zoom into the bound entangled area shown in (a). 




FIG. 2: (Color online) Same as Fig. [T] for nonvanishing values of 7 in Eq. (9|. Also in these cases a perfect agreement is 
observed between the areas detected by the lower bound and the exact areas defined by Eqs. (llOll. 



the figures), the corresponding expressions read 

/Sent = -17a +197 



/3i 



PPT 



±3 v^l + 6 a - 39 a2 + 30 7 - 102 a 7 + 33 7 
1 

16 

±3\/4- 12 a- 15a2 - 12 7 



^(-2+lla+ll7 



66 a 7 — 157-^ 

(10) 



and are indicated respectively by the dashed eUipses and 
full line in Figs, [l] and |2] the gap region between the 
ellipses and the full line, inside of the PPT region, defines 
the bound entangled area. 

Given the perfect agreement between the results ob- 
tained using the quasi pure approximation and those 
from optimal witnesses, we now address a class of states 
which hitherto could not be characterized by the latter. 
These are states that extend "beyond lines" in the above 
sense, i.e., which cannot be generated by application of 
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FIG. 3: Parameter space of the "beyond line" states as de- 
fined in Eq. (Ill for 7 — —1/12. Tlie color code is as before. 



only one Weyl operator. We choose the following family: 

p = (1 - a - /3 - 7)i + 7P00 + aPio + I3Pqi . (U) 

We can once again easily scan the whole region of param- 
eters for a fixed 7, and new areas of bound entanglement 
are found, as illustrated in Fig. [3] 

Our above results suggest that bound entangled states 
in general occupy a finite volume in the even higher 
dimensional state space of all 3 x 3 states, as was ac- 
tually proven in |21j . and underpinned for random bi- 
partite states of dimension 2 x 4 in [32]. In a similar 
vein as in |23j . we now show numerical evidence that 
our approach to detect bound entangled states is ro- 
bust, i.e., it does not only work for Bell-diagonal states, 
but also for a finite bound entanglement volume around 
them. As an example, we choose a "line-state" (pbc) with 
{a,/3,7} = {-0.092,0.04,0.2148} (and Cqp = 0.018) as 
defined in ([9]) above, and mix it with a Hilbert- Schmidt- 
distributed random mixed state (phs)i as follows: 



P 



= (l-e)pbc + epHS, withO<e<l. (12) 



In this way we are able to explicitly construct a ball of 
bound entangled states. 

As before, the bound entangled fraction is identified by 
intersection of the area with positive partial transpose 
and that with non-vanishing quasi pure approximation 
Cqp. The result is illustrated in Fig. [4| by the distribu- 
tion of Cqp and of the negativity Af [6P(in the inset) , as 
function of the variable e. Similar plots are obtained for 
different initial bound entangled states, even for "beyond 
line" states. Note that Af vanishes if and only if the state 
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FIG. 4: Distribution of quasi pure concurrence Cqp, Eq. 
and of the negativity A/" (inset) of 3 x 3 random states, gen- 
erated by admixture of Hilbert-Schmidt-distributed random 
mixed states pns to the bound entangled reference "line" 
state pbe, according to Eq. (|12[). Saturation indicates the bin- 
populaton of the histogram for a fixed e. 



has positive partial transpose, and can thus be used to 
demarcate the associated parameter range. 

For each value of e, 1000 states pns were randomly 
chosen, such that for e = 1 we recover the entanglement 
characteristics of our sample - where both, Cqp and M, 
exhibit a broad distribution. In the opposite limit, e = 0, 
Af vanishes identically, while Cqp — 0.018 (identifed by 
the dashed line) - indicating that pbc is bound entangled. 
The figure shows that already for e < 0.01 all sampled 
states have a positive partial transpose (vanishing Af) 
but non- vanishing Cqp, and are thus bound entangled. 
Since the states around Pbe are randomly chosen and are 
greater in number than the dimension of state space, the 
probability of all of them lying in a hyperplane is zero 
and thus they explore all directions in state space - the 
convex hull of these points forms a body of finite volume 
in state space. Therefore, our numerical result explicitly, 
and effortlessly, spots a finite volume of bound entangle- 
ment in the state space. Nevertheless, we cannot rule out 
the existence of separable states in the constructed vol- 
ume. But given the large number of our sample, and the 
fact that separable states form a convex set, the proba- 
bility of such event is vanishingly small. 
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